Improving CMB non-Gaussianity estimators using tracers of local structure 
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Local non-Gaussianity causes correlations between large scale perturbation modes and the small 
scale power. The large-scale CMB signal has contributions from the integrated Sachs Wolfe (ISW) 
effect, which does not correlate with the small scale power. If this ISW contribution can be removed, 
the sensitivity to local non-Gaussianity is improved. Gravitational lensing and galaxy counts can 
be used to trace the ISW contribution; in particular we show that the CMB lensing potential is 
highly correlated with the ISW signal. We construct a nearly-optimal estimator for the local non- 
Gaussianity parameter /nl and investigate to what extent we can use this to decrease the variance 
on /nl. We show that the variance can be decreased by up to 20% at Planck sensitivity using galaxy 
counts. CMB lensing is a good bias-independent ISW tracer for future more sensitive observations, 
though the fractional decrease in variance is small if good polarization data is also available. 



I. INTRODUCTION 

Non-Gaussianity is a possible signature of early- 
universe physics which should be measured to high ac- 
curacy by precision CMB observations [l|. Forthcoming 
data from the Planck satellite can constrain the bispcc- 
trum with local shape to ±5 (see below), which may be 
comparable to the expected signal if inflation produces 
purely Gaussian fluctuations [2j. Any improvement in 
the error bar is very welcome, since the current limits 
of /nl = 32 ± 21 [l[ are already constraining /nl to be 
comparable to the Planck error bar. One promising ad- 
ditional constraint comes from scale-dependent bias in 
large-scale structure surveys @. In this paper we con- 
sider a different possibility: we investigate to what extent 
the estimates of local non-Gaussianity from CMB data 
alone can be improved by using a tracer of the large-scale 
matter distribution. 

Local non-Gaussianity results in a non-zero bispectrum 
of specific shape: it causes correlations between the large- 
scale fluctuations and the small-scale power. For exam- 
ple, a large cold spot on the CMB may have more fluc- 
tuations on small scales than over a large hot spot. Non- 
Gaussianity can be produced by several effects: it could 
be present in the primordial fluctuations when they are 
generated during inflation, being a powerful discrimina- 
tor of different inflation models; it will be generated by 
evolution of the perturbations between generation during 
inflation and last-scattering; and it can also be generated 
by gravitational lensing between the last-scattering sur- 
face and our observations. The latter two effects are ex- 
pected to be present, and in principle can be accurately 
predicted and subtracted off the observed signal to re- 
cover constraints on the primordial contribution. 

The precision of non-Gaussianity constraints is lim- 
ited by observational noise and resolution. Generally the 



larger the number of small-scale modes that can be ob- 
served, the better the constraint, but this is limited from 
Planck at I ~ 1600 where the noise becomes important. 
On smaller scales, secondaries may also be an important 
source of confusion. In this paper we seek to improve the 
constraint not by improving the number of small-scale 
modes, but by increasing the signal in the correlation of 
the small-scale modes with the large scale modes. 

The large-scale CMB temperature perturbation has 
contributions both from last-scattering — which are ex- 
pected to correlate with the small-scale modes at last- 
scattering if there is local non-Gaussianity — but also 
from the integrated Sachs Wolfe (ISW) effect 0. The 
ISW effect arises from red- and blue-shifting of CMB pho- 
tons as they move through evolving potentials along the 
line of sight, with the induced temperature anisotropy 
given in terms of the Weyl potential \& by the line-of- 
sight integral 
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where the dot denotes a conformal time derivative. The 
contribution is mainly from redshifts z < 3 when dark 
energy starts to effect the growth of structure, and hence 
probes fluctuations a long way from the last-scattering 
surface at z ~ 1000. The ISW contribution to the large- 
scale CMB is therefore expected to be uncorrelated to 
the small-scale signal at last scattering, even if there 
is local non-Gaussianity. The ISW contribution effec- 
tively acts as a source of noise on any /nl estimator 
using only the observed CMB temperature. If this con- 
tribution to the temperature could be subtracted off we 
would be able to infer the large-scale temperature at last- 
scattering, which would then be better correlated with 
the small-scale power, giving better constraints on lo- 
cal non-Gaussianity because of the absence of the ISW 
'noise'. A similar idea has recently been applied to sta- 
tistical anisotropies using large-scale structure data 
Here we also consider using information in the lensing- 
induced CMB trispectrum to reconstruct the lensing po- 
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tcntial, and then use this to subtract the ISW contribu- 
tion from the observed temperature. 

The paper is organized as follows. In Section |n] we 
briefly review the key theoretical concepts behind the 
construction of estimators for /nl- In Section IIIII we 
derive the estimators we will be using in this paper, 
whose aim is to remove the effects of the ISW on es- 
timators for /nl- In the first sections of the paper we 
only consider CMB data, comparing idealized models of 
the Planck satellite and a possible future CMB mission 
(specifically EPIC: Experimental Probe of Inflationary 
Cosmology 0). In Section ITVl we discuss the instrumen- 
tal noise and reconstruction noise if using a CMB lcnsing 
reconstruction. Fisher forecasts for the improvement in 
error bar using estimators are presented in Section [Vj 
along with analysis. In Section IVIII we briefly consider 
how to use multiple tracers, and finally in Section IVIIII 
we present our conclusions. 

Throughout the paper we adopt a ACDM cosmology 
with parameters £7 m = 0.25, SI a = 0.75 and Ho = 73km 
s" 1 Mpc -1 . 



II. THE CMB BISPECTRUM 

The bispectrum, the three-point function, is a useful 
statistic for the detection of non-Gaussianity. The CMB 
bispectrum -B/iZ 2 / 3 is defined for a statistically isotropic 
universe by 
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where the quantity in parentheses is the Wigner 3j- 
symbol and the ai m are the spherical-multipole coeffi- 
cients of the CMB. It will be non-zero due to CMB lcns- 
ing and other non-linear effects, but is most studied as 
a probe of primordial non-Gaussianity. If there are mul- 
tiple observed fields, {q^}, for example the CMB tem- 
perature and polarization, or some probe of the matter 
density, we can also define the more general bispectra 

/Si) Sj) Jk) \ _ B ijk ( h h h \ ,*s 
\ ( ii imi ( ii 2m 2 q hm 3 / ~ r> hhl 3 \ T n 1 m 2 m 3 J [ ' 

= b?£ la J dQY hmi Y hm2 Y l3m3 .(4) 

In the second line we defined the reduced bispectra t , 
which can only be defined if l± + 12 + 13 is even, as it must 
be for scalar and gradient (E-mode) fields. 

The amplitude of a CMB bispectrum of known shape is 
denoted /nl, and we shall focus on bispectra due to local 
primordial non-Gaussianity where most of the signal is 
in 'squeezed triangles' (the correlation of one large-scale 
mode with two small-scale modes). The amplitude of 
any non-Gaussianity is already constrained to be small so 
we can use an Edgcworth expansion about the Gaussian 
distribution in order to motivate optimal estimators Q. 



This allows a weakly non-Gaussian full-sky PDF to be 
expressed as a sum of its cumulants: 
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where in this expression the q/ m represents a vector of 
fields, for example just the CMB temperature, or a com- 
bination of several fields e.g. q; TO = (ai m ,ipi m ) where 
ipi m is a matter tracer field that is correlated with the 
ISW. The matrix Cj = (q/ m q|, n ) is the correspond- 
ing covariance. We can perform the differentiation in 
Eq. ([5]) , and after applying selection rules for the Wigner 
3j-symbol and discarding monopole terms (see Q for full 
details) we obtain a simplified version of the PDF, 
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Maximizing the likelihood, rfl "f = 0, gives an estima- 
tor for the bispectrum amplitude Jn l of the form 



h<h<h 



where each bispectrum term is estimated using 
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and we have introduced a permutation factor Ai ± i 2 i 3 
which is I when l\ ^ I2 ^ h, 6 when l\ = I2 = h and 
2 otherwise. In this paper we will use the Fisher 'ma- 
trix', F, to quantify the error <t/ nl = F _1//2 on /nl (the 
amplitude of B^[ ol ), which is given by Q 
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For the case of a single field, the full sky Fisher error 
is determined simply by 



6 .^r! Ci x Ci 2 Ci 3 
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Since the estimator is derived in the limit of small non- 
Gaussianity, the Fisher error estimate is also only valid 
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in this limit; if significant non-Gaussianity is detected 
the non-Gaussian variance can significantly modify the 
result [|[ . In this paper we focus on expected error limits 
in the null-hypothesis that there is no non-Gaussianity. 



III. NON-GAUSSIANITY ESTIMATORS 

In this section we derive and analyse different ways 
to implement estimators for /nl- We will examine esti- 
mators that not only include information from the CMB 
temperature a; m , but also from the a tracer of the po- 
tential field responsible for the ISW (ipim), and E-mode 
polarization information (labelled 'E'). 

A. Incorporating a tracer of the ISW 

The simplest possible way of removing the ISW is a 
'subtraction estimator', in which an estimate of the ISW 
contribution to the CMB temperature is simply sub- 
tracted from the observed map. This is essentially the 
same procedure as used in Ref. [U when trying to study 
CMB anomalies at last scattering, using observed galaxy 
number densities as a tracer for the ISW. We define an 
ISW-cleaned temperature anisotropy as, 
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where tj)i m are the multipole coefficients of some tracer 
field and Cf = (ai m ipi m )i with analogous definitions 
for Cf^' and Cj T . Here we have assumed that we want 
to subtract all CMB temperature that is correlated with 
the tracer; this may not be quite correct since very large- 
scale perturbations anti-correlate the last-scattering and 
ISW signals at the 10% level on large scales. 

If we use qi m = ai m in Eq. (J9j) , then we will obtain 
an error for the subtraction estimator. The result will 
be the same as Eq. (fTt)|) . except we must replace C; with 
Q = (ai m a* m ), given by 

(cf*) 3 
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In addition -Bj 1 i 2 z 3 must also be changed to reflect 
the fact that we are now interested in the bispectrum 
(ai imi ai 2m2 ai 3m3 ) rather than (ai imi ai 2m2 ai sm3 ) . 

The subtraction estimator is a suboptimal way of com- 
bining a measurement of tpi m with the CMB tempera- 
ture. The optimal way to include the extra information 
is to use the vector of fields q; TO = (ai m ,ipi m ) T ', where 
the expected error is determined by Eq. © and we in- 
clude the information from all eight possible bispectra 
(TTT, TTip, TtpT, tjjTT, Ttpip, tpTtp, 4>i)T, t/jipi/)) and the 
covariance is 



C, 



(13) 



Alternatively, instead of using the vector of fields q; m = 
(ai m ,ipim) T 'j we could do an equivalent optimal anal- 
ysis using a pair of orthogonalized variables qj m = 

{ai m ,^i m ) T where (ai m ipi n ) = 0. To the extent that 
the tracer and the ISW-cleaned temperature probe in- 
dependent parts of the universe (at very different red- 
shifts), we expect the de-correlated fields di m and ipi m 
to be independent as well as uncorrelated. In this ap- 
proximation there are only two non-zero bispectra (B aaa 
and B^^), and hence the bispectrum estimated from 
the ISW-cleaned temperature is only suboptimal to the 
extent that it is neglecting information contained in the 
Qiptpip bispectrum. If we only wish to use ipi m on large 
scales, this extra information should be small since there 
are only a small number of modes. Also due to compli- 
cated non-Gaussian properties of any likely tracer, it may 
also be a good idea not to include this less reliable infor- 
mation, in which case using the subtraction estimator is 
nearly the best thing one can do. We will later provide 
a quantitative comparison. 



B. Including polarization Information 

Using an analogous method to the inclusion of infor- 
mation from an ISW tracer, it is also possible to include 
the effects of polarization by using a two-component vec- 
tor q; m = (ai m , Ei m ) T where Ei m are the multipoles of 
the E-mode polarization 0, [lfj. This optimal T-E es- 
timator is the same as the optimal T-ip estimator, but 
replacing if) with E im . 

We can further combine CMB polarization and tem- 
perature information with a tracer i\) of the ISW. There 
are several ways to include E and ip information: the 
simplest is to compute an estimator which we will label 
TE - this is the same as the TE estimator except in- 
stead of using CMB anisotropics ai m we use ISW-cleaned 
anisotropics a; m . The other more optimal alternative is 
to construct the optimal T — E — ip estimator. To do this 
we can use the vector of fields q; TO = (ai m ,ipi m , Ei m ) T , 

where the terms (C -1 )*"' in Eq. ^ are now individual 
terms taken from the 3x3 covariance, 
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Alternatively we could consider a subtraction estimator, 
using the vector of fields (T, E) where 



Ei, 



Ei, 



(15) 



which is by construction uncorrelated to tp. Note that 
although the ISW does not contribute significantly to the 
E-polarization, there is nonetheless a correlation due to 
a correlation between the large-scale E-mode signal from 
reionization with local (z ~ 2) structures (see Fig. [1] and 
Ref. [HI for details). 
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FIG. 1: The correlation between the CMB lensing potential ip 
and ISW (solid line), and the correlation between the lensing 
potential and the CMB temperature (dashed) and polariza- 
tion (dash-dotted) for a standard ACDM model with reioniza- 
tion optical depth r = 0.09. The lensing potential and ISW 
are very well correlated. 



IV. CMB LENSING AND NOISE 

A promising tracer of the ISW is the CMB lensing 
potential ipi m , which is given in terms of the line of sight 
Weyl potential 'F by 



tp(h) 



dx 
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where Vn'0 gives the deflection angle, 770 — x is the confor- 
mal time at which the photon was at position %n, fx (x) 
is the comoving angular-diameter distance, and the CMB 
is well approximated by a single source plane at comoving 
distance %*. In concordance ACDM models the lensing 
potential coincidcntally happens to have a very similar 
kernel to the ISW, as indicated by the ~ 90% correlation 
between the lensing potential and the ISW as shown in 
Fig. [1] The lensing potential can be reconstructed us- 
ing the statistical anisotropy induced in the small-scale 
CMB by lensing (see Refs. QJ, QJ| for a review). 
Error bars on /nl depend on the noise properties of the 
measuring instrument - we will consider Planck and a 
possible EPIC configuration in this paper. Instrumental 
noise will be important both in determining the accu- 
racy of our lensing reconstruction, and in determining 
the errors on the measurement of the temperature and 
polarization anisotropies. 

We approximate the instrumental noise as isotropic, so 
it contributes a term to the power spectrum iVj, with 
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TABLE I: Relevant noise parameters for Planck and EPIC for 
a frequency of 150 GHz. a 2 characterizes the detector noise 
in dimensionless units (Eq. (|17l> ) and Afwhm is the beam full- 
width half-maximum. 





o 2 


#fwhm /arcmin 


Planck 


2.69 x 10" 17 


7.1 


EPIC 


8.65 x 10~ 21 


5.6 



where a 2 is the white detector noise power and 

$fwhm is the beam full-width half-maximum 
(FWHM). In this paper we focus on simple models of 
the Planck and EPIC experiments, with parameters 
summarized in Table HI For EPIC we have used the 
information provided in Q for the 150 GHz frequency 
band, assuming other frequencies are used only for 
foreground subtraction. For simplicity we shall assume 
a full-sky observation with isotropic noise, since our 
purpose in this paper is to assess whether ISW cleaning 
is potentially useful rather than describing a realistic 
analysis. 

Using galaxy counts to trace the ISW, good data 
should yield a ip which is essentially cosmic variance lim- 
ited on large scales. However, when using CMB lensing 
reconstruction ip to trace the ISW, there may be signif- 
icant reconstruction noise. If we reconstruct the lensing 
potential using temperature information alone, following 
[l5j and using the same notation, this noise can be mod- 
eled by adding a noise term to Cf given by 



A L = (2L + 1) 
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where Cf f = Cj T + N t and, 
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+ Sj0] , 



and Ei = I 2 + I. 

We can also reconstruct the lensing potential using in- 
formation from polarization, which we will henceforth 
refer to as an E-B reconstruction. In this case, the noise 
term to be added is, 



A L = (2L + 1) 
where in this case, 
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For low experimental noise, the E-B reconstruction is 
much better than the temperature-only reconstruction 
since there are expected to be no unlcnscd B-modes on 
small scales, so the observed small-scale B modes are 
probing lensing directly without confusion from an un- 
lcnscd signal. The quadratic estimator reconstructions 
considered here are somewhat suboptimal once the re- 
construction noise becomes small [16| |. and iterative or 
optimal estimators can do significantly better. 



V. TRACING THE ISW SIGNAL USING CMB 
LENSING ALONE: RESULTS AND ANALYSIS 

We use CAMB [TtJ to calculate the required transfer 
functions, correlation matrices and local bispectra. The 
results for the estimators introduced in Section IIIII are 
presented in Table |TT] for both Planck and EPIC noise 
parameters. 

To quantify the possible theoretical improvement in 
/nl error bars using ISW subtraction we can consider 
the hypothetical case in which we know, a priori, the ex- 
act form of the ISW, and hence can subtract it perfectly 
from the CMB maps. For the temperature-only estima- 
tor this scenario is referred to as 'T (ISW = 0)' in Table 
HI! showing that the error bars could potentially be im- 
proved by about 10% (~ 20% decrease in variance). This 
is not dramatic, but nonetheless can be considered sig- 
nificant if compared with the cost of observing longer to 
correspondingly reduce the small-scale noise. Since the 
lensing potential is highly correlated with the ISW, we 
would expect a perfect reconstruction of the lensing po- 
tential ip to be close to this ideal result, and the results 
in Table [TT] show that indeed the improvement remains 
at nearly 10% if the lensing potential could be measured 
perfectly on large scales. 

We can compare these ideal cases with more realistic 
possibilities shown in Table [TT1 corresponding to the esti- 
mators we discussed in the preceding sections with noise 
on the lensing reconstruction. Results for all estimators 
including the ip field are quoted with the inclusion of a cut 
in I - all ipim terms with I > 50 are discarded so that only 
large scales are being included, rather than also includ- 
ing non-Gaussian signals intrinsic to ip that in practice 
are likely to be untrustworthy due to the complicated 
statistics of ip (a full joint analysis of the primordial and 
lensing bispectrum and trispectrum is beyond the scope 
of this paper, but could potentially improve constraints 
further). The results are insensitive to the precise value 
at which we take the I cut. 

Where relevant we calculate the error on /nl both with 
and without noise from the lensing reconstruction. We 
reconstruct the lensing potential using both temperature 
information and polarization information. 

In the case of Planck, the noise in the lensing recon- 
struction is sufficiently large that there is little improve- 
ment in the /nl error using the lensing potential as the 
ISW tracer. For EPIC the reconstruction is much better: 



TABLE II: Errors on /nl derived using various estimators 
for both Planck and EPIC noise parameters. The estima- 
tor labels are explained in the text. Where relevant we have 
indicated whether or not we have included noise in the re- 
construction of the lensing potential. We have also identified 
how we have calculated the reconstruction noise - whether it 
be from temperature (T) or polarization (E-B). The results 
for the T tp estimator are quoted for a cut in / excluding terms 
at I > 50 



Estimator 


Error on /nl 
Planck 


EPIC 


T 


5.90 


4.74 


T (ISW = 0) 


5.32 


4.28 


Tip (no reconstruction noise) 


5.39 


4.31 


Tip (T reconstruction noise) 


5.80 


4.60 


Tip (E-B reconstruction noise) 


5.86 


4.35 


Subtraction (no reconstruction noise) 


5.41 


4.34 


Subtraction (T reconstruction noise) 


5.86 


4.69 


Subtraction (E-B reconstruction noise) 


5.86 


4.39 


TE 


5.19 


2.44 


TE (no reconstruction noise) 


4.92 


2.36 


TE (E-B reconstruction noise) 


5.19 


2.39 


TE (T reconstruction noise) 


5.19 


2.42 


TEf/> (no reconstruction noise) 


4.90 


2.35 



including noise on the lensing reconstruction we can still 
reduce the noise on /nl estimated from the temperature 
by ~ 8%. We could also use temperature and polariza- 
tion data to reconstruct the potential - this will result in 



a further small improvement. Using results from [14| we 
can estimate that using temperature together with po- 
larization the reconstruction noise could be reduced by 
a further factor of 10, leading to a noise on /nl of 4.32 
(an improvement of a further 0.03 from the case where 
we only use polarization information). 

Excluding lensing reconstruction noise we see that the 
subtraction estimator is close to optimal (the error on 
/nl is close to that of the Tip estimator). The subtrac- 
tion estimator has the added advantage of being fast to 
compute, much faster than the 'optimal' estimators, and 
thus in practice may be preferable. 

Using polarization information significantly improves 
the error on /nl (by a factor of about two [1, Hoj ) 
even without ISW subtraction. At EPIC sensitivity 
more small-scale temperature and polarization modes 
are available, so the relative importance of the largest- 
scale modes is lower than for Planck. However, even for 
Planck, polarization can in principle help significantly, 
since it probes somewhat different triangles because of 
the phase shift between the polarization and temperature 
transfer functions, and also provides another handle on 
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the large-scale modes. In both cases, if polarization in- 
formation is used the fractional improvement from using 
ISW subtraction becomes smaller. In reality the large- 
scale polarisation data may be hard to determine due to 
sky cuts and foregrounds; if we exclude polarisation data 
at I < 20 the TE estimator will perform 3% better than 
the TE estimator, so ISW subtraction still gives some 
improvement. As mentioned in Section Hill we could also 
consider a subtraction estimator using the vector fields 

(f, t) , which improves the error by only a further < 1% 

compared to the more basic TE estimator. 

The best possible estimator we can construct is one 
that optimally includes information from temperature, 
polarisation and the ISW effect: this is the TE-0 estima- 
tor in TablcQI] We see that in the ideal case for EPIC this 
optimal combination reduces the error by about ~ 4% 
compared to using TE alone, and with Planck the im- 
provement is at the ~ 6% level. The TE subtraction 
estimator is another (simpler) way of combining temper- 
ature, polarisation and ISW information - as we can see 
from Table HH these subtraction estimators are almost 
optimal, and have a comparable error to the TE^ esti- 
mator. 

As an additional test of the results we set the ipi m to 
equal the exact value of the ISW (rather than tracing 
it by the lensing potential). As expected the results in 
Table [TT| change by < 1% - this is further confirmation 
that the lensing potential is an extremely good tracer of 
the ISW. 



VI. SCALE-DEPENDENCE 

The ISW effect is a large-scale phenomenon, and thus 
is more important at low I. Thus we would also expect 
that cleaning the ISW signal from the temperature im- 
proves the signal to noise mostly in triangles with one 
very low-l side. If we want to test different sources of 
a detected non-Gaussianity signal, resolving any scale- 
dependence will be very useful pil UM i an d ISW-cleaning 
may then be relatively more useful to improve the con- 
straint on triangles with the lowest- 1 on one side. 

We analyse the ratio of two errors, the optimal T-only 
estimator and the Tifj estimator, restricting triangles to 
have at least one side lower than a certain I threshold. For 
clarity we excluded lensing reconstruction noise. As can 
be seen from Fig. [21 as the I threshold is reduced the frac- 
tional improvement becomes larger, reflecting the fact 
that ISW cleaning is most useful for large-scale modes. 
We performed the same analysis for estimators including 
polarisation, with the same conclusion - the ratio of the 
errors for the TE estimator and TE estimator are plotted 
in Fig. 1 

In reality, we only wish to subtract ISW contributions 
to T that are generated by local structures. The correla- 
tion matrix C T ^ however includes the contribution from 
large-scale modes, which stretch from last-scattering to 



z < 3 which may contain additional information on non- 
Gaussianity. So for example, in the Tip subtraction es- 
timator, we ought to use C T ' sw -<°™ 1 ' 4 ' instead of C T ^' . 
Doing this marginally improves the result of the subtrac- 
tion estimator, but given that the subtraction estimator 
has already been shown to be close to optimal, this im- 
provement is < 1%. 

VII. OTHER TRACERS OF THE ISW - 
GALAXY NUMBER COUNTS 

As discussed in Section HVl the lensing potential is a 
very good tracer of the ISW. However, we have seen that 
the noise on lensing reconstructions is high at Planck sen- 
sitivity, so other tracers may be much more useful until 
future CMB missions such as EPIC are able to recon- 
struct a cosmic-variance limited large-scale lensing po- 
tential. 

We focus here on galaxy number counts. It is well 
known that measurements of galaxy densities can be 
used to probe the matter density on large-scales, and 
can therefore be used as a tracer of the ISW effect (see 
e.g. Ref. [13, [H| and references therein). Compared 
with using gravitational lensing as a tracer, the situa- 
tion with galaxies is more complicated since the galaxy 
density is generally a biased tracer of the matter distribu- 
tion and hence potentials. For standard Gaussian mod- 
els, the bias is expected to be nearly scale-independent 
on large scales at a given redshift, and hence should make 
a reliable tracer for the ISW. In the presence of primor- 
dial non-Gaussianity the situation is more complicated 
however, since the modulation of the small-scale power 
spectrum by large scale modes gives rise to strongly scale- 
dependent bias on large scales [3j]. A full joint analysis of 
scale-dependent bias is beyond the scope of this paper; 
instead we will assess the use of the large-scale galaxy 
distribution as a tracer of ISW under the null hypothe- 
sis that the primordial non-Gaussianity is negligible, so 
that the bias is scale independent. As with CMB lensing, 
here we are only interested in the large-scale part of the 
matter density as a probe of ISW; using the small-scale 
galaxy density as a direct probe of non-Gaussianity is 
potentially promising but very much more complicated 
and currently unproven. 

We assume galaxy number counts trace the syn- 
chronous gauge matter density perturbation 5 and use 
them to reconstruct the ISW potential, using a model 
for the galaxy redshift distribution (see @) given by, 

^ <xz a atp\~(z/z*f] • ( 22 ) 
dz L J 

We take a = 2, /3 = 1.5 and we vary z*. We assume full 
sky coverage and that shot noise is unimportant on the 
scales of interest (low I). Taking the source distribution 
given by Eq. (|22[) , we can calculate the new Cf ^ and C ; T * 
using CAMB and re-evaluate the subtraction estimator. 
The results are shown in Table IIIII 
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FIG. 2: To investigate scale-dependence we calculate the Tip estimator and the T-only estimator (left-hand figure) and the 
TE and TE estimators (right-hand figure), restricting the sum to include triangles with one side lower than the I threshold 
plotted on the x-axis. The y-axis is the ratio of the errors as calculated using the two estimators (Tip/T and TE/TE). The solid 
line represents Planck and the dashed line corresponds to EPIC. We see that cleaning the ISW signal from the temperature 
improves the signal to noise mostly in triangles with one very low-Z side, as the ISW is predominantly a large-scale phenomenon. 



TABLE III: Subtraction estimator errors for Planck and EPIC 
when using galaxy number counts to trace the ISW effect. 

z, Planck EPIC 

0.2 5.63 4.51 
0.7 5.47 4.37 
1.5 5.69 4.55 



For Planck, number counts with sources going out to 
z ~ 1 are likely to be the best ISW tracer since the 
lensing reconstruction noise is large. For EPIC, CMB 
lensing reconstruction performs comparably, and has the 
advantage of directly probing the same potentials rather 
than having to model bias. Galaxy number counts will 
perform better when there is a close match between the 
form of Eq. (|22|) and the ISW kernel - as we see from Table 
IIIII the error estimates can vary significantly depending 
on the parameters used in Eq. (|2"2"j) . Current low redshift 
data such as 2MASS offers almost no improvement. 

We can easily combine information from multiple trac- 
ers or redshift bins. Labeling two tracers as ipi and ip2, 
we may first evaluate the ISW-cleaned CMB anisotropics 
using ipi, ai m , as before, 

ai m - a hn ^ . (23) 

°; 

Then we clean ai m again using the field ip2 to give our 
final result a\ m , 

aim = aim - ^ 2 ■ (24) 



In general, if we have information from TV fields ipi we 
can construct a vector, v = {Tisw, ipi, ip2, V'iv) an d 
define a covariance C = (vv T ). Then the maximum like- 
lihood estimator for the cleaned temperature given the 
information from the ipi is (indexing vectors and arrays 
from zero) is, 

aim = ai m + }^ /^-n ■ ( 25 ) 

i=l ^ >oo 

This defines an optimal linear combination of the differ- 
ent probes which gives the best tracer of the ISW. 

VIII. CONCLUSIONS 

The principal conclusions from this work are: 

1. We used the Edgcworth Expansion to derive op- 
timal estimators that take into account combina- 
tions of the CMB temperature, an ISW tracer, and 
CMB polarization. We showed that a simple easily- 
computed 'subtraction estimator' can be used to re- 
move an estimate of the ISW contribution to large 
scales, and that this is close to optimal. 

2. If the ISW could be removed perfectly, the 
maximum reduction in the error on local non- 
Gaussianity /nl is at the 10% level. In practice 
any realistic tracer of the ISW will do worse than 
this. 

3. The CMB lensing potential is an excellent tracer of 
the ISW, with a correlation close to unity over the 
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entire range of I applicable. In the noise-free case, 
we found that we can effectively remove the effects 
of the ISW, approaching the optimal limit. For 
Planck the lensing reconstruction is too noisy to 
significantly improve the /nl error, however with 
a possible future satellite (EPIC) we can still re- 
move the effects of the ISW at a near-optimal level 
- the reduction in the error on temperature-only 
/nl estimators is ~ 9%, and around 5% when the 
polarization information is included. 

4. The ISW is a scale-dependent effect, so ISW- 
subtraction is most useful for constraining triangles 
with one very low-Z side. It may be most helpful 
for testing scale-dependent non-Gaussianity mod- 
els, where the CMB provides the only information 
on the largest scales. 

5. For Planck, using galaxy counts to trace the ISW 



signal is much better than using CMB lensing re- 
construction; we estimate a ~ 9% reduction in the 
error on /jvl if the ISW kernel is well matched. For 
EPIC, CMB lensing would work well and provide a 
robust alternative to using galaxy tracers, though 
the improvement is small if large-scale polarization 
can be used reliably. 
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